The quantum ground state of self-organized atomic crystals in optical resonators 
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Cold atoms, driven by a laser and simultaneously coupled to the quantum field of an optical 
resonator, may self-organize in periodic structures. These structures are supported by the optical 
lattice, which emerges from the laser light they scatter into the cavity mode, and form when the 
laser intensity exceeds a threshold value. We study theoretically the quantum ground state of these 
structures above the pump threshold of self-organization, by mapping the atomic dynamics of the 
self-organized crystal to a Bose-Hubbard model. We find that the quantum ground state of the 
self-organized structure can be the one of a Mott-insulator or a superfluid, depending on the pump 
strength of the driving laser. For very large pump strengths, where the intracavity intensity is 
maximum and one would expect a Mott-insulator state, we find intervals of parameters where the 
system is superfluid. These states could be realized in existing experimental setups. 

PACS numbers: 



I. INTRODUCTION 

Self-organization in systems of atoms and light has 
been observed in numerous experiments, some of which 
are reported in Refs. [H-@|- Several experimental real- 
izations show the formation of spatially-ordered atomic 
structures, which organize in the potential they form [l|, 

i| . The basic mechanism behind the observed dynam- 
ics can be summarized by considering that the refrac- 
tive index of the atomic medium, which is related to the 
atomic spatial density, is itself determined by the light 
fields via the mechanical effects of atom-photon interac- 
tions. Hence, the effective dynamics the atomic center 
of mass undergoes is determined by potentials and/or 
forces, which in return depend on the atom position and 
velocity distributions Q- 

In this context, one remarkable example is the forma- 
tion of regular patterns of atoms in the standing wave 
of a high-finesse optical cavity, arising when the atoms 
are transversally driven by lasers. Experimental signa- 
tures are the phase-locking between the field at the cav- 
ity output and the driving laser phase, according to two 
possible values with difference equal to 7r. This occurs 
when the laser frequency is smaller than the cavity-mode 
frequency and its intensity exceeds a threshold value @ . 
The phenomenon, first predicted by numerical simula- 
tions according to a semiclassical model for atoms and 
light [8|, can be understood in terms of the atoms be- 
ing trapped in the potential which originates from the 
light coherently scattered from the laser into the cavity 
mode. The largest intracavity amplitude is supported by 
the atomic configurations in which all atoms scatter in 
phase into the resonator, corresponding to the patterns 
where the interatomic distance is a multiple of the wave- 
length of the cavity field. Since the possible configura- 
tions which fulfill this condition in a standing- wave cavity 



are two, one shifted with respect to the other by half a 
wavelength, the difference between the possible phases of 
the emitted field is exactly tt. A further theoretical study 
discussed the phenomenon in terms of second-order phase 
transition and determined the pump threshold within a 
semiclassical, mean- field model @. 

Additional novel phenomena arise in this system when 
the quantum mechanical properties of light and matter 
are relevant. Theoretical works studied the dynamics 
of entanglement between atoms and fields during self- 
organization in a similar setup [Tol - fT3 | . In Ref. [l3| the 
pump threshold for self-organization at T = was deter- 
mined when the atoms are assumed to be forming a Bose- 
Einstein condensate. Most recently, the onset of self- 
organization has been experimentally observed in a Bosc- 
Einstein condensate coupled to an optical resonator 
In such situation, an open question regards the nature of 
the quantum state of the system when the pump inten- 
sity is increased well above the self-organization threshold 
value. In this regime, in fact, the height of the poten- 
tial increases, and one would expect localization of the 
atoms at the potential minima when the atom number is 
increased and hence the field in the cavity is larger. On 
the other hand, interparticle interactions due to s-wave 
collisions compete with such dynamics, favoring situa- 
tions in which the number of atoms per site is small. 
This framework is reminiscent of the paradigmatic Mott- 
insulator / superfluid quantum phase transition [lij , and 
in particular its realization with ultracold atoms in op- 
tical lattices In the case here considered, however, 
the dependence of the potential on the atomic distribu- 
tion (and in particular, the expectation that its depth in- 
creases with the atomic density) makes it a priori unclear 
whether incompressible states can exist at large laser in- 
tensities. 

In this article we study the quantum ground state of 
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this self-organized system, in the setup sketched in Fig.Q] 
In particular, we determine under which conditions in- 
compressible states exist at large pump intensities. For 
this purpose, we extend a procedure, developed for an 
atom-cavity system in which the cavity mode is pumped 
by a laser [16l - ll8l |. and derive a Bose- Hubbard model 
for the atoms trapped in the potential their generate, 
whose coefficients depend on the atomic density. Using 
the strong-coupling expansion [19j . we identify the pa- 
rameter region where incompressible states of the self- 
organized atomic pattern may exist. This study finds 
particular motivation from experimental progress, which 
achieved the coupling of ultracold atoms with the optical 
mode of high-finesse cavities jf| [20l - [25l | . Indeed, we will 
argue that incompressible self-organized states could be 
observed in the setups of existing experiments [f| [26| . 
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FIG. 1: (color online) Setup of the system, showing matter- 
wave density distribution (blue region) in the effective poten- 
tial created inside an optical resonator by coherent scattering 
of laser photons (red line). Here, is the Rabi frequency, giv- 
ing the strength of the coupling between the atoms and the 
laser field, A is the wavelength of laser and cavity mode, and re 
is the rate at which photons leak from the cavity mode. The 
quantum state of the system could be inferred by measuring 
the light at the cavity output, or by Bragg spectroscopy using 
a weak probe (see Sec. |V)| , 



This article is organized as follows. The equations at 
the basis of the physical model are introduced in Sec. lU 
and the relevant approximations leading to the deriva- 
tion of the effective atomic dynamics in the light poten- 
tial are extensively discussed. In Sec. Mil we assume 
that the atoms are tightly bound at the minima of the 
cavity potential inside the cavity, and derive an effective 
Bose-Hubbard model, whose coefficients are analytically 
determined within a modified Gaussian ansatz. The pa- 
rameter regimes, where incompressible states of the self- 
organized atomic gas are found, are derived in Sec. IIVI 
within the strong coupling expansion method [19]. Con- 
clusions and outlook are presented in Sec. |V] while the 
appendix provides further details complementing the dis- 
cussion in Sees. Mil and IIVI 



II. THE MODEL 

Identical bosonic atoms, with mass m, are confined in- 
side an optical resonator, which is assumed to be a high- 
finesse cavity sustaining well resolved modes. The atoms 
are prepared at ultralow temperature T and interact with 
each other by means of s-wave scattering. Furthermore, 
the atoms couple to light via a dipolar transition between 
the electronic ground and excited state \g) and |e) at the 
optical frequency ujq. In particular, the atomic dipole is 
simultaneously driven by a laser and by the field of one 
cavity mode, both linearly polarized and whose wave vec- 
tors are perpendicular to one another, according to the 
geometry of the setup sketched in Fig. [TJ We assume the 
atoms to be strongly confined in the plane perpendicu- 
lar to the cavity axis, here chosen to be the x-axis, so 
that the motion is essentially along x. In absence of the 
resonator, the atomic dynamics for the atoms in states 
\g) and |e) are governed by the Hamiltonian "Ho + T~L egi 
where 



(1) 



9 j (x) 



with gj the strength of the state-dependent collisional 
interaction, and where H e g describes s-wave scattering 
between atoms in different atomic states. The Hamil- 
tonian in Eq. (p} is written in second quantization, 
where operators *f?i (x) and "jf (x) annihilate and cre- 
ate an atom at the position x and in the electronic state 
i = g,e. The atomic field operators obey the bosonic 

commutation relations 9i(x), 9'Jx') = Sij5 (x — x') 

and [9i(x),9j(x')] = 0. 

The cavity is assumed to be a high-finesse resonator 
sustaining well-resolved modes, of which the mode at fre- 
quency uj c is (quasi) resonant with the atomic transition. 
We denote by a and the annihilation and creation 
operators of a cavity photon at frequency uj c . The coher- 
ent coupling between atoms and light is described by a 
Hamiltonian in the Rotating- Wave approximation, which 
reads 



hA a / dx^l (x) * e (x) - aA c a f a 



H dx (g cos (kx) a + Q(x)) *| (x) * g (x) + H.c. 



and which is here reported in the reference frame ro- 
tating at the laser frequency u> p . Here, Q(x) denotes 
the strength of the coupling to the laser field (Rabi fre- 
quency), go is the cavity- mode vacuum Rabi frequency, 
cos(fcx) gives the spatial dependence of the cavity-mode 
function with wave vector k, while A a = uj p — uio and 
A c = ujp — lj c denote the detuning of the laser frequency 
from the atom and from the cavity-mode frequencies, re- 
spectively. The coherent dynamics of coupled atoms and 
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electromagnetic field are then governed by the full Hamil- 
tonian 

T~L = Ho + T~L\ . 

This far we just considered Hamiltonian dynamics. 
We now introduce the sources of noise and dissipation, 
which in the present model are spontaneous radiative de- 
cay of the the electronic excited state |e) at the lifetime 
7 and cavity losses at rate 2k due to the finite trans- 
mittivity of the mirrors. Noise is here described within 
the Heisenberg-Langevin formalism, and the Heisenberg- 
Langevin equations for the atomic and field operators 
read QJ 



(3) 



-i (.go cos (kx) a f + Q(x)) * e (x) - ^fl(t)V e (x) 

* e {x) = ~ [* e (a:) , Ho] + (iA„ - |) * e (x) , (4) 
-i^ g (x) (go cos (kx) a + fi(x)) + /y* s (x) f- m (t) , 

+V2Ka in (t). (5) 

Here, the noise operators fi n (t) and a- m (t) are the 
Langevin forces, with vanishing mean value and 
<//„(*)/*(*')> = (4n(tM')) = 0, while (fUWLW) = 
(a in (t)al n (t')) = S(t - t'), see for instance Ref. [271 ]. 

A. Effective dynamics 

We now assume a time-scale separation, such that the 
electronic variables relax to a local steady state on a much 
shorter time scale than the one characterizing the dy- 
namics of external and cavity degrees of freedom. This 
regime for the cavity degrees of freedom requires that 
the parameters satisfy |A + ij/2\ 3> \Cl\,goVN, |A C |, k. 
Moreover, a change in the external degrees of freedom 
can be neglected over the typical time scale of the internal 
degrees of freedom provided that |A a + i-y/2j kbT/Tl, 
with kb Boltzmann's constant. In addition, we assume 
|A„| 3> \fl(x)\, goy/N 7/2, so that we can identify 
a time scale in which coherent evolution is taking place 
while dissipation due to spontaneous decay can be dis- 
carded. In these limits Eq. (j4]) can be set to zero and at 
leading order gives 

V e (x) ~ * s (x) [ago cos (kx) + fl(x)\ / A a , (6) 

such that Eq. ([5]) takes the form 

a — ~kcl + i [A c — Uoy] a — iSqZ + V2~Kai n , (7) 

where Uq = 5o/Aa is the maximal depth of the single- 
photon potential and So — ,go^/A a is the maximal am- 
plitude of scattering a laser photon into the cavity mode 



by a single atom. In writing Eq. ([7J we have assumed, 
moreover, that Sl(x) — Qf(x), with |/(x)| < 1. The spa- 
tial distribution of atoms and fields is now contained in 
the operators y and Z, which read 

y = J dxcos 2 (kx)^l(x)^ g (x), (8) 

Z = J dxf(x)cos(kx)^l(x)^ g (x) , (9) 

and act on the Hilbert space of the atoms. The opera- 
tors y and Z are Hermitian and commute, [y, Z] = 0, 
as they both depend solely on the atomic density. They 
are the quantum analogous of the semiclassical "bunch- 
ing parameter" B — J2j cos,2 (k x j)/N and of the "spa- 
tial order parameter" 8 = J2j cos(kXj)/N, respectively, 
which are defined for an ensemble of N atoms at the po- 
sitions X\, . . . ,1.2V) see for example [9j. These parameters 
characterize a semiclassical mean-field description of self- 
organization of the atoms in the cavity field, such that 
for O — >• ±1 the atoms are in a self-organized state at 
the even or odd antinodes of the cavity standing wave, 
while B provides their degree of localization about these 
points [28| . 

The cavity field can be eliminated from the matter- 
wave equations assuming a second characteristic time 
scale determined by the rate |A C + ik\, in which the 
cavity field approaches a local steady state. This re- 
quires that the effective coupling strength between mat- 
ter waves and cavity field is smaller than this rate, namely 
|A c +i«| 3> \Sq\Vn, \Uq\N, and that the uncoupled mat- 
ter wave dynamics is slower, |A C + hc| 3> kqT '/%. In this 
limit, one finds a ~ a^°\ with 



,(°) = 



SoZ 



2Ka ln (t) 



(a c - u y) + iK 



(10) 



where we have used that the operators in the numerator 
commute with the operators in the denominator. Using 
in Eq. ([3]) the solutions for f e (i) and a, Eq. ([5]) and 
Eq. (|I0[) . at lowest order in perturbation theory we find 

i> g (x) = — ^7 [V g (x) , Ho] - iU cos 2 (IfcrJo^V' ' 
-iSof(x) cos(kx)a m -f g - iS f(x) cos(kx)t> 3 a (0) 

n 2 



-i—f(xy* g 



(11) 



Here is now a function of the atomic density and 
hence does not commute with the field operator *f> g (x). 

B. Mechanical effects of the cavity field 

Equation (TTOl) for the cavity-field operator gives 
the field state at leading order in a perturbative expan- 
sion, in which it is assumed that the cavity field follows 
adiabatically the matter-wave dynamics. Furthermore, it 
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is seen that it is a function of the atomic density, which 
enters both in the numerator and in the denominator of 
the expression. In this work we will be interested in deter- 
mining the atomic quantum ground state, which emerges 
from a coupled dynamics between cavity field and atoms, 
giving rise to the mechanical potential sustaining such 
state. 

In order to gain further insight into the effect of the 
cavity potential on the atoms, we consider the semiclas- 
sical limit, which can be found from Eq. (TTT1) . In this 
limit, we identify the terms which explicitly depend on 
the atomic position with mechanical forces. The corre- 
sponding potential takes the form 

so 2 f( T \2 

V(x) = V 1 cos 2 (kx) + V 2 cos(kx) + AW_,(12) 

with Vi = aC7 <a (0)+ a (0) ) and V 2 = 2hS Re{(a^)}. The 
third term on the right-hand-side of Eq. (|12[) is the a.c- 
Stark shift due to the laser field. In the following we 
assume that the external laser field drives the atoms uni- 
formly and the Rabi frequency fl does not depend on the 
position x (i.e. f(x) — 1). The Fourier transform of the 
potential is hence composed by a term at wavevector 2k 
whose amplitude is proportional to the number of pho- 
tons, and which is due to the dispersive coupling between 
photons and atoms. This term corresponds to a potential 
with periodicity A/2, analogous to an optical lattice in 
free space. Differing from an optical lattice in free space, 
the number of photons depends on the atomic density 
via the " bunching parameter" . The second component 
is at wave vector k and thus oscillates with periodicity 
A. It originates from coherent Raman scattering of laser 
photons into the cavity field and its amplitude is hence 
proportional to the cavity electric field amplitude, but it 
also depends on the positions of the scatterers through 
the "spatial order parameter" 0. The two Fourier com- 
ponents give rise to an effective potential with spatial 
periodicity A as displayed in Fig. [T] We remind that 
potential (|12[) is found in a semiclassical limit. In this 
limit, the expression we find agrees with the ones derived 
in [H, Hl| within a semiclassical model. 

The model of a semiclassical potential for the atoms 
does not take into account dissipative effects, which may 
be due for instance to cavity decay (while dissipation due 
to back-action of the resonator over the atoms is here a 
priori neglected as we assume the adiabatic limit). In- 
deed, in the dynamics of the atoms we find non-hermitian 
terms, which are characterized by the position-dependent 
coefficient 

j'(x) =lm{S Q a (0) }cos(kx), 

scaling with the linewidth k of the resonator. Corre- 
spondingly, input noise gives rise to fluctuations in the 
height of the potential, which are described by the op- 
erator eti n (i) in Eq. (TT0|) . The effect of cavity quantum 
noise will be here discarded, assuming that the detuning 
I A c | 3> k (more precisely, |A C | ^> k, \ Uoy\, which is con- 
sistent with the assumption of adiabatic elimination of 



the cavity degrees of freedom). In this regime, the atoms 
move in a dispersive potential, which is produced by the 
dynamical Stark coupling with the cavity mode photons 
and by the Raman-scattered field of the laser. 

C. Tight-binding limit 

The motion and steady state of the semiclassical model 
discussed in the previous section has been studied in sev- 
eral theoretical works. Numerical simulations demon- 
strated that the system self-organizes, such that the 
atoms localize at the antinodes of the cavity field, where 
their coupling is maximum, according to an array with 
periodicity A This behavior appears when the laser 
intensity exceeds a critical value [8j, |9| , and has been con- 
firmed experimentally Q . Localization of the atoms may 
take place around two possible set of minima, one cen- 
tered at the positions x 2 °j — jX (where cos(kx^j) = +1), 
with j integer, and the other at the positions = 

(2j + l)A/2 (where cos(kx^j +1 ) = — 1), while the cav- 
ity field is a coherent state with amplitude proportional 
to the pump [29l . |30| . In the limit in which the atoms 
are well localized at the minima one can expand poten- 
tial (fT2")l till second order in the fluctuations of the parti- 
cles position. For a particle localized at xf^ the potential 
takes the form 

V(x) ~ V w + V {2) x 2 (13) 
where x — Xj — Xj and = V(Xj), while 

V '"> = + {U " {2) + ( - 1>,(A < - 

(14) 

where the potential due to the other particles enters in a 
mean-field approach through the mean value of operators 
Z and y. Note that when the atoms are self-organized at 
the set of positions {x^ '}} (respectively, {a4j+i})i tnen 
the sign of (Z) is positive (respectively, negative). 

Consistently with the approximations made so far, the 
dynamics can be characterized by the dispersive dipolar 
potential provided that | U (y) | , | U (Z) | < |A C |, which 
allows for the adiabatic elimination of the cavity field 
from the matter- wave equations. This implies that, for 
the regime we consider, the sign of this term of the po- 
tential is solely determined by the sign of the detuning 
between cavity and pump. One hence finds that the po- 
sition x^ are minima provided that A c < 0. This be- 
havior can be understood in terms of positive feedback of 
the system, which is warranted whenever the sign of (Z) 
is opposite to the one of the harmonic potential V^ 2 \ 
and hence when A c < 0. This latter condition leads 
to the property that the conservative force, due to the 
potential, is attractive [3l|. It is a sufficient condition, 
provided that the detuning A Q between atom and laser 
is negative, A a < (and consequently Uq is negative), 
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as it is visible from the numerator of Eq. (fl"4)) (note that 
for the situations here considered, |(Z)| > (y), see also 
Ref. [Hj]). Indeed, in this regime the minima of the po- 
tential, formed by light scattering of the atoms, are also 
minima of the standing-wave potential of the cavity field, 
so that the attractive forces of both contributions add up 
to confine the atoms. On the contrary, when A a > the 
opposite situation is realized: the minima of the poten- 
tial due to light scattering are now maxima of the cavity 
standing-wave potential: the term proportional to Uq in 
the numerator of Eq. (fT4"f has opposite sign than the 
detuning A c , and the overall effect is to make the ef- 
fective potential, resulting from the two contributions, 
shallower. An extensive discussion on the semiclassical 
forces as a function of the detunings can be found in 

Refs. mmi. 

Finally, the amplitude of the semiclassical detuning de- 
pends nonlinearly on the atomic density through quantity 
A c — U y in the denominator. Bistability effects due to 
this nonlinearity have been discussed for the situation, 
in which the cavity is pumped by a laser [T3, [H, [25|, |H, 
l32l l33j . In this article we will focus on the regime where 
I A c — Uo(y)\ 3> k, when the potential is conservative 
and dissipative effects can be neglected. This parameter 
regime is far away from the bistability region. Corre- 
spondingly, this term will be treated as a small correc- 
tion of the semiclassical potential in the framework of 
perturbation theory. 



III. EFFECTIVE DYNAMICS IN THE CAVITY 
POTENTIAL 

A. Derivation of the effective Hamiltonian 

We now assume that the atoms are tightly confined 
in one set of minima of potential (fT2"j). say, at the posi- 
tion X2j = jA. This assumption is clearly based on the 
semiclassical dynamics, as a potential cannot be simply 
singled out from Eq. Our aim is to identify the 

quantum ground state by considering the full quantum 
dynamics. For this purpose, starting from the assump- 
tion of tight confinement, we perform a Wannier decom- 
position of the atomic field operator, 



*9 ( X ) =^2 W j( X ) b 3 > 



(15) 



where operator bj annihilates a particle at the site cen- 



tered in x ^j j and Wj(x) is the Wannier function, taken 



to 



be real valued and with center in x^j ■ This decomposi- 
tion is based on the assumption that the atoms are in the 
lowest band of the semiclassical potential, which is justi- 
fied at ultralow temperatures. Within the decomposition 
of Eq. (HU, Hamiltonian {TJ takes the form Ho ~ Hq h , 
with 

H* H = E Q N + E X B + ~uJ2m (m - 1) - Mo^V , (16) 



where N = bjbj is the atom number operator, and 
the coefficients read 



B = J2j (bjbj+i + H.c.^ is the hopping operator, while 
fficients read 

f ( h 2 d 2 \ 

Ei = dxw t ( x )l- — — Wi+i ( X ) (17) 



U 



gin / dxwi(x) 4 



(18) 



and g\D = g g is the scattering strength in one dimen- 
sion. In writing Eq. (|16p we made the nearest-neighbor 
approximation, and introduced the chemical potential fi, 
assuming a grand-canonical ensemble. We remark that 
the Bose-Hubbard expansion of Hamiltonian ([T5]) relies 
on localization of the atoms at the minima of the poten- 
tial inside the cavity. 

We now use Eq. (fT5|) in Eq. {JTTJ) , and determine the 
equations of motion for the operators bi by multiplying 
with wi (x) and integrating over the position. We find 

bi = -\[buH$ H ] -iU Q A^[J k + Ji^M (0) 
- iSo^W [Z b t + Z x Bi] - iS [Zok + ZiBi] -4 (0) 

.n 2 



(19) 



which is written in the regime in which |A C | 3> |E/o!V|- 
Here, B t = b l+1 + b^ = [b h B], and A (0) ~ a' ', such 
that 



^ (0) = SojZoN+^B) ;1 



-^L^(J N + .hB)) , 

A c + IK / 

which considers only nearest-neighbor coupling and 
where we used 



Z 

y 



ZoN + ZiB, 
J N + J X B, 



with 



Z( = J dx Wi(x) cos (kx) Wi + i(x) , 
J i = I dx Wi(x) cos 2 (kx) Wi+e(x) . 



(21) 
(22) 



(23) 
(24) 



and £ = 0, 1. At first order in the expansion in the small 
parameter {/ /|A c |, Eq. (fl"9|) can be exactly cast in the 
form 



(25) 



with 



ttBH 

CQED 



ns 2 



{ZoN + ZiB) 



U 



Z 



(26) 



x (Z J N 3 + Zo.hNBN + J Q Z X (N 2 B + BN 2 ) + o(Z 2 ))] , 
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where we consider k <C |A c | and we truncated the second 
term in an expansion at first order in \ZiUqN/A c \ and 
| JiIIqN/ A c \. Note that in expression (f26f the operators 
are symmetrically ordered, so to obtain Eq. (ITT))) . 

The total dynamics of the system is now rendered by 
the Hamiltonian 



off 



H, 



BH 



H, 



BH 
CQED 



(27) 



where the first term of the right-hand-side, defined in 
Eq. ([T5)l . describes the hopping between sites, due to the 
quantum fluctuations, and the on-site interaction emerg- 
ing from s-wave scattering, while the second term con- 
tains the contributions due to the mechanical effects of 
the laser and cavity potential, which determine and sus- 
tain the atomic pattern. We note that Hamiltonian (|26|) 
is similar to the Bose-Hubbard Hamiltonian, however, 
its terms are nonlinear in the number and hopping op- 
erators. Moreover, it differs also from the usual Bose- 
Hubbard model due to the appearance of a term propor- 
tional to the squared of the hopping operator, B 2 . 

The nonlinearity in Eq. (|26p is due to the coupling with 
the cavity field, which scales with the number of atoms 
and which determines the confining potential: the larger 
is the atomic number the stronger is the coupling. This 
effect is exquisitely due to cavity quantum electrodynam- 
ics, and originates from the constructive interference with 
which the atoms scatter the laser photons into the cavity 
mode. We also note that, when considering higher order 
terms in Uo/ A c , an exact effective Hamiltonian -which is 
a function solely of operators N and B fulfills Eq. (fH))) - 
cannot be found unless one resorts to a suitable thermo- 
dynamic limit (l8j . In fact, these terms emerge from the 
nonlinear coupling between photons and atoms inside the 
resonator, in which the atoms act as a refractive index. 
The effect is hence highly nonlocal, and it is significant 
when the strong coupling regime is warranted. An imme- 
diate consequence is that Wannier functions will depend 
on the particle number. Most remarkable is the fact that 
the Hamiltonian is quadratic in the hopping operator, 
even when the nonlinearity deriving from the atomic re- 
fractive index, C/o3^, is negligible. This property is due to 
the collective scattering of the atoms into the resonator. 
Nevertheless, it does not affect substantially the prop- 
erty of the ground state at large potential depth (small 
tunneling rates) as we will show. 



B. Mapping to a Bose-Hubbard model 

For the purpose of studying the quantum phases of 
the many-body system, we define a thermodynamic limit. 
This is here so defined, that the density of the atoms is 
kept fixed as the particle number N and the cavity-mode 
volume go to infinity. Denoting by 2K the total number 
of sites of the lattice, such that in the considered con- 
figuration there are K lattice cells, the cavity coupling 
strength is inversely proportional to the squared root of 
the cavity mode volume, and thus scales as go oc l/y/K. 



The other parameters in Eq. ([27|) scale as Uq — uq/K 
and So = sq/\/K, where uo and so are constants, while 
the number operator is now N — Khq, with no the den- 
sity operator, giving the number of atoms per site. This 
scaling is such that, as K — > oo, the parameters in Hamil- 
tonian (|27l) depend only on the atomic density. We also 
note that in this thermodynamic limit a Bose-Hubbard 
form can be also derived when the terms scaling with Uq 
are taken into account, see for instance [l8[ for a discus- 
sion. In this thermodynamic limit we rewrite Hamilto- 
nian (|2"T|) in the Bose-Hubbard-like form, 



H = -t 1 [no]B+t 2 [n Q ]B 



where 



2 , U[no] 



>J ni (rtj - 1) - /i[n ] N , 
(28) 



2Z Z 1 + n ^{Zlh + 2J Z Q Z 1 ) 



*a[no] = [1 + O(n /A c )] , 



»o 



M n o] = Ho-E - n Q h—Z 1 + no — J 



and where U, Zj, and Jj are functions of no through the 
Wannier functions. Insight into the dependence of the 
coefficients on the density can be gained using a (modi- 
fied) Gaussian ansatz in place of the Wannier functions, 
as reported in Appendix Q] 

Differing from the ordinary Bose-Hubbard model [Tij . 
in the case of self-organized atomic patterns the coeffi- 
cients entering the effective BH Hamiltonian do depend 
on the atomic density. This property implies that the 
atomic densit y is not solely determined by the chemical 
potential fl7llT8j. We also observe that, in addition to 
the nearest-neighbor hopping arising from the kinetic en- 
ergy, there is a hopping term due to the cavity potential 
which contains a term proportional to B 2 . This latter 
contribution describes long-range correlations among the 
atoms mediated by the scattered photons. Close inspec- 
tion shows that, in the tight-binding regime, its coeffi- 
cient is of higher order with respect to the coefficient 
multiplying B, see App. Q] Within the considered ther- 
modynamic limit, this term gives rise to a small correc- 
tion to the spectrum of the lowest excitations and corre- 
sponding eigenstates of the Mott-insulator states, which 
must, however, be taken into account when determining 
the ground state of the system for larger values of the co- 
efficient Z\ , scaling the tunneling coefficient t 2 and some 
terms of t\. 



IV. INCOMPRESSIBLE PHASE OF 
SELF-ORGANIZED ATOMIC PATTERNS 

We now study whether the self-organized atomic pat- 
tern exhibits an incompressible ground state. Therefore, 



(29) 
(30) 

(31) 
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we focus onto the regime in which self-organization has 
set in, namely, when the laser amplitude exceeds the 
threshold value [H, HH and the laser frequency is smaller 
than the cavity-mode frequency, A c < 0. The atoms are 
assumed to be localized at the set of minima {x 2 j = jA}, 
such that Zq > 0. 

In order to display the phase diagram for this system, 
we consider the chemical potential and the amplitude of 
the pump field sq, which indirectly determines the tun- 
neling coefficients. In fact, for sq — > oo the intracavity 
field potential is deepest and tightly localizes the atoms 
at their minima, and the tunneling parameters t\ and t 2l 
Eqs. (|29|) and (|30)l . vanish correspondingly (as one can 
check by considering the dependence of the coefficients 
Zj and Jj on so, see Appendix 0]). 

W e ap ply the strong coupling expansion developed in 
Ref. [191 ] in order to determine the size of the incompress- 
ible states as a function of the pump amplitude and of the 
chemical potential. For this purpose we start by consid- 
ering the limit sq — > oo, assuming the system to be in a 
Mott-insulator state with an integer occupation per site 
uq. We calculate the free energy of the Mott-insulator 
state Em (no) and of the defect states E±(no), which are 
obtained by adding or removing a particle to the Mott- 
insulator state. The energy of the defect state is found 
by means of degenerate perturbation theory, expanding 
in power of the small parameter ti/U up to third order. 
In the parameter regime we consider, this expansion cor- 
responds also to a first-order expansion in the parameter 
t 2 . The phase boundaries between the Mott-insulator 
and the compressible phase correspond to the situation 
where the energy difference between the defect state and 
the Mott-insulator state vanishes, namely, when 



E ± (n o )-E M (no) = 0. 



(32) 



The values of the chemical potentials /t*i ) and fJ>r no y 
at which the Mott-insulator phase with integer occupa- 
tion riQ becomes unstable when adding and removing a 
particle, respectively, fulfill the equations 



t 2 

L l 2 
— Tin 

u 



E + (n ) - E M (n ) = -M(„ ) + Un - 2ii(n + 1) + 
t 3 

+ j^n (no + l)(n + 2) 
+ t 2 2(2n + l)(n + 2) =0, (33) 
E-(n ) - E M (n ) = (J,^ - U (n - 1) - 2txn 



t 2 t 3 
+ — (n + l) 2 + - 

+ t 2 2(2no + l)(no-l)=0. (34) 



Mno + l) 2 + 77U0K + l)(no - 1) 
U (7 Z 



When fi+ no) 



> 



then uf 
• 



and a, determine 



the upper and lower boundaries of the parameter region 
where the system is in the Mott-insulator state. The 
phase diagrams reported afterwards are computed using 
the analytical formulae of the strong coupling expansion, 
Eqs. (|33p and (|34[) in the thermodynamic limit, fixing the 
values of the parameters uq and sq. The coefficients are 



evaluated using a modified Gaussian ansatz in place of 
the Wannier functions, see Appendix for detail. 

Figure [2] displays the phase diagram for the considered 
system, showing three incompressible, Mott-insulator 
phase regions with n = 1,2,3, for A c = —50k and 
uo = —0.1k. When plotted as a function of l/s , the 
Mott-insulator regions become larger as a function of 
the laser amplitude as the number of atoms increases. 
This result can be understood with simple arguments, 
in fact the trapping potential is constituted of the light 
coherently scattered by the atoms: As the number of 
atoms increases, the cavity-field intensity, for a given 
pump amplitude, is larger, thus the potential becomes 
deeper and the configuration more stable. Another re- 
markable property is the appearance of gaps between the 
Mott-insulator lobes. In the limit of vanishing tunneling 
coefficients t\,2 — > (corresponding to so — > 00), this 
behavior implies that the degeneracy among states with 
different densities, found in the ordinary BH model, is 
here removed. This behaviour results from the depen- 
dence of the single-particle, on-site energy on the aver- 
age density. In particular, while inside the lobes the den- 
sity per site no takes integer values, in the gap between 
the lobes the density changes continuously in the interval 
r < no < r + 1 (with r = 1, 2, 3, . . . integer number), as 
shown in Fig. [3] The region outside the Mott-insulator 
lobes is presumably a peculiar superfluid state, which en- 
tails photonic excitations. A study of its properties close 
to the self-organization threshold has been presented in 
Refs. The appearance of the gap is understood, 

when one considers the dependence of the coefficients on 
the density. In fact, at so — > 00 (corresponding to putting 
the tunneling elements to zero), the Mott-insulator re- 
gion with no — 1 is found for values of fx such that 
< Hi < U%, where ui — u[no — 1] and U\ — U[no — 1] 
(see Eq. (fT5|> ). Similarly, the region at no = 2 is found 
for U 2 < u 2 < 2U 2 (with U 2 = U[n Q = 2]). The rela- 
tion U\ < U 2 (Uj < Uj + i) leads to the appearance of 
parameter regions where the system is compressible even 
at zero tunneling. This relation is a consequence of the 
fact that, for a fixed amplitude of the laser field, the in- 
tracavity field becomes larger as the number of atoms 
(scatterers) is increased. 

Figure @] displays different phase diagrams, obtained 
for various values and signs of the detuning A a , keeping 
A c fixed. We first focus on the case in which A a < 0, 
when the potential due to the dynamical a. c. -Stark shift, 
induced by the coupling with the cavity mode, adds up 
to the potential due to light scattering, giving rise to a 
stronger attractive dipole force. In this case, we observe 
that for larger values of uq the gap between the chemical 
potentials of Mott-insulator states at different densities 
increases. Moreover, the larger is uo, the larger is the 
area covered by Mott-insulator states at larger densities: 
incompressible states are found at smaller values of the 
laser amplitude. The opposite tendency is found when 
A a > 0, and hence when the potential due to the a.c- 
Stark shift competes with the scattering potential. In 
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this case, as |uo| is increased, the gap between the Mott- 
insulator lobes becomes smaller, while the area covered 
by Mott-insulator states with larger densities decreases. 
In particular, incompressible states are observed at larger 
values of the laser pump amplitude. This behavior can be 
understood from the dependence of the tunneling coeffi- 
cient, Eq. (1291) . on the various terms: The term propor- 
tional to Mo scales with the onsite density, and it hence 
becomes more important as no is increased. For values 
of the detuning, such that the cavity potential adds up 
to the scattering potential (A a < 0), at a fixed laser 
intensity the confinement is larger for a larger density 
of atoms. Correspondingly, the parameter region where 
the system is compressible becomes larger. On the con- 
trary, for A a > the cavity potential tends to make the 
scattering potential shallower, the confinement tends to 
depend less strongly on the density and the gap between 
the Mott-insulator lobes is smaller. 



1 1 1 1 




n = 3 






n o = 2 


\ 






\ 1, : 




1 ' 





FIG. 2: (color online) Phase diagram, showing the Mott- 
insulator states at different densities no = 1, 2, 3 in the plane 
reporting the chemical potential fi — fl/Ui, Eq. (|3ip . and the 
inverse of the laser amplitude, 1/sq (in units of cavity loss rate 
ft). The parameters are A c = —50k and uo = —0.1k. The 
arrow indicates the threshold value of the laser amplitude, for 
which the system self-organizes [l3| ]. The chemical potential 
is here reported in units of the on-site energy Ui , correspond- 
ing to the value of Eq. (118[) for the density of one atom per 
site, with <7id/(-EbA) = 3.74 • 10 -4 and Er the recoil energy, 
considering a gas of 87 Rb atoms whose dipole-transition, at 
wavelength A = 830 nm, is coupled to the cavity mode. 



We remark that, while in this model atom- atom inter- 
action is essentially s-wave scattering, it should be con- 
sidered that, when the atomic occupation of each site 
exceeds unity, then dipole-dipole interaction, including 
superfluorescence, will be relevant. In our model we ne- 
glect collective effects in the spontaneous decay outside 
of the cavity. In particular, for the parameter regime 
we consider, incoherent processes are suppressed for the 
large atomic detunings we choose. On the other hand, 
one can identify parameter regimes where coherent cou- 
pling between the atoms at a site may be relevant, and 
which should be included in an effective onsite interaction 
term. 

We finally estimate the parameter regime required for 
observing our predictions, considering the setup reported 




0.5 - 

q| 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 3 3.5 

FIG. 3: (color online) On-site density no as a function of the 
chemical potential ji (in units of the on-site energy U\) for 
so — > oo (corresponding to ti,tz — > 0). Same parameters as 
in Fig. |U 
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K/S 

FIG. 4: (color online) Mott-insulator states at no = 1, 2 in the 
plane Ji/Ui-k/sq for A c = —50k and (i) uo = —2k (red solid 
line), (ii) uo = 0.1k (green dashed line), and (iii) uo — 10k 
(blue dotted line). The other parameters are the same as in 
Fig.0 



in Ref. [26[. Here, ultracold 87 Rb atoms are loaded in- 
side a resonator, where a mode couples quasi-resonantly 
with the transition at wave length A = 780 nm. The 
parameters are go/ix = 14. 1 MHz, k/2-k = 1.3 MHz and 
7 /2tt = 3 MHz. Taking |A a |/27r = 60 GHz, fl/27r = 60 
MHz, we find U /2n ~ 0.03 MHz and S /2tt = 0.1 MHz, 
while the effective rate of spontaneous emission is of the 
order of few KHz. For a number of atoms N ~ 10 4 , 
the relation SqV~N ^> K is fulfilled, warranting the exis- 
tence of a time-scale in which the atoms may experience 
a conservative potential formed by their scattered pho- 
tons. The relation |A C | ^> \Uq\N,k is satisfied taking 
|A c |/27r = 100 MHz. For these parameters, the value 
So/ft = 10 at which incompressible phases are found in 
Figs. [2] and HI corresponds to S VN ~ 10k. This latter 
value corresponds to the one found for the experimental 
parameters in [26| , showing that the regime, where quan- 
tum effects in self-organized atomic patterns are visible, 
can be accessed with existing experimental setups. 
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CONCLUSIONS 



relevant to the dynamics of atom-photon interactions. 



The quantum ground state of a self-organized atomic 
pattern, confined by the potential of the field scattered 
into the cavity mode, can be incompressible provided 
that the laser intensity exceeds a certain critical value, 
which is found above the threshold value at which the 
system self-organizes. In addition, intervals of values of 
the chemical potential are identified, in which the state 
is compressible even when the tunneling vanishes. These 
appears as gaps between the Mott-insulator lobes in a 
phase diagram we derive, and their origin can be under- 
stood from the competition between confinement, which 
is steeper at large atomic density for a given laser inten- 
sity, and interparticle collisions, whose effective strength 
also depends on the density via the localization of the 
atomic wave function in the intra-cavity field poten- 
tial. Our analysis is based on an effective Bose-Hubbard 
model, according to the assumption that the atoms are 
tightly bound at the minima of the intra-cavity field po- 
tential. Since the height of such potential depends on the 
atomic density, the coefficients of the model depend on 
the atomic density giving rise to the peculiarities encoun- 
tered in the phase diagram. Comparison with typical pa- 
rameters of experimental setups, dealing with ultracold 
atoms inside of optical cavities, shows that the prediction 
here made could be experimentally tested in existing se- 
tups. 

Detection of the quantum properties of the system 
may be performed by measuring the light using a probe 
field [f| l34l . IHI , in a direction of emission corresponding 
to a Bragg scattering angle. A spectral analysis of the 
probe intensity will provide information on the state of 
the system, which can be studied as a function of the laser 
pump intensity, and hence of the tunneling rate [3(1 |37| . 
Sweeping the phase diagram in the direction jl may be 
performed in various ways. One possibility is to imple- 
ment a shift of the atomic ground state by means of a 
magnetic field or by off-resonant coupling with a third 
electronic state. 

In this work we focused on the quantum ground state of 
the self-organized structure, arising by driving the atoms 
with a strong, transversal laser field. The mechanical 
potential results from the scattering potential, with peri- 
odicity equal to the wavelength A, and from the potential 
due to the cavity field mode, with periodicity equal to the 
wavelength A/ 2. The regime we considered is the one in 
which the scattering potential is the strongest. In this 
regime, it is expected that the light at the cavity out- 
put is coherent, while the intensity is independent on the 
number of atoms when the strong-coupling regime is war- 
ranted [29[, 30]. An interesting question is what are the 
properties of the emitted light when the cavity-field mode 
potential dominates. In this regime, in a semiclassical 
model it has been found that the light at the cavity out- 
put can be squeezed, and in general anti-bunched [38l , l39| . 
An intriguing question is whether and how these proper- 
ties are modified, when the atomic quantum statistics is 
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1. Coefficients in the Gaussian Approximation 

The solution of the Bose-Hubbard model requires the 
evaluation of the coefficients, which are integrals con- 
taining Wannier functions. The Wannier functions are 
solutions of the Schrodinger equation for a single parti- 
cle in presence of the potential in Eq. (|12l) , which itself 
depends on the spatial distribution of the atoms inside 
the cavity. The problem is hence nonlinear. We solve 
it by combining an analytical and a numerical approach. 
In first instance, we replace the Wannier functions by 
modified Gaussian functions, Wi(x) — » Wi(x), with 

*iO) = (l + 9i(x) - |(1 + \P) {g i+ i(x) +gi-i(x)) 



(35) 



and 



(x-x;)^ 



(36) 



is a normalized Gaussian function of width a. The 
form of the modified Gaussian function in Eq. (|35[) war- 
rants orthonormality up to second order in <5, where 
S = J dxgi(x)gi + i (x) = e _7r "^ fc "°' ) and where the over- 
lap between second neighbors is of the order 0(<5 4 ). We 
remark that the coefficient in (|35p at second order in 
5 warrants the correct normalization. The width a of 
the Gaussian function g+ [x) is found by minimizing the 
single-particle energy E, for the Schrodinger equation in 
presence of potential (IT21 . For later convenience we in- 
troduce the dimensionless parameter 

V = k 2 a 2 , (37) 

which is y <C 1, consistently with our assumptions. In 
terms of the parameter y this consists in finding the min- 
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imum of the equation 



E(y) = ^ + F(y)-^ly., 



(38) 



where Er — H 2 k 2 /2m is the recoil energy and we intro- 
duced the functions 



F(y) = Vi + V 2 , 
G(y) = V 1 +V 2 /2, 



where V\ and V 2 are 
% 



hU S 2 Z(y) 2 
K 2 + (A c - U Q y(y)) 2 ' 
2hSjZ(y)(A c -U Q y(y)) 
At 2 + (A c - U y(y)) 2 



(39) 
(40) 



(41) 
(42) 



Note that the cos-functions in Eq. (fT2j) have been ex- 
panded up to second order in the displacement from the 
minima. Minimization with respect to y leads to the 
nonlinear equation 



y 2 = 



Er 



2F'{y) - yG'{y) - G{y) ' 



(43) 



with F'(y) = dF(y)/dy and G'{y) = dG(y)/dy, and 
which is solved numerically by recursion. 



Using the value which minimizes the single-particle en- 
ergy, y := j/min, we evaluate the coefficients we need in 
order to study the problem. They are reported below as 
a function of y: 



E 

Jo 
Zo 
Ei 



Er 



1 + S Z 



An* 



2y 
i + c-y 

2 ' 
(1 + A5 2 )e-y/ 4 



y 



Ji = o, 



25e-y' i [1 + -5 2 
4 



U = gi D k 



1 + 3<5 2 - 4ji 
\/2vry 



(44) 

(45) 
(46) 
(47) 
(48) 
(49) 

(50) 



where Er = h 2 k 2 /2m is the recoil energy and we dis- 
carded terms of higher order than S 2 . Note that the 
coefficients Zj are here reported for the pattern localized 

at the positions x^j = jX. 
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